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The Transient Response of a Dipole Antenna Using a 
Summation Approach in Technique (SAT) 


Mohamed N. Morsy, Senior Member, IEEE and Walter K. Kahn, Life Fellow, IEEE 


Abstract—A novel summation approach in technique (SAT) im- 
proving on the usual computation by multiple reflections to find 
the transient response was presented in a previous paper. This tech- 
nique leads to accurate answers consisting of relatively few terms. 
The transient response at large values of time can be found without 
need to know explicitly and laboriously add all multiple reflections 
for the whole preceding range of time involved in producing the 
response. In this paper, the transient response for a realistically 
modeled pulse generator exciting a dipole antenna modeled as a 
distributed parameter load is investigated using a generalization 
of the SAT. A detailed comparison between this technique and the 
standard superposition of multiple reflections is presented. 


Index Terms—Computational electromagnetics, EMC, light- 
ning, pulsed radiation, time-domain techniques, transients. 


I. INTRODUCTION 


T is well known that time- domain techniques are very 

useful to compute the response of structures excited by 
transient electromagnetic signals. The transient response to 
generated pulses for different kinds of loads has been studied 
extensively in the literature [2]-[12]. The transient response 
may be formulated, in general, as an infinite number of suc- 
cessively delayed waves. It may be computationally difficult to 
accurately predict the transient response at larger values of time 
by the conventional method of adding all multiple reflections 
involved in producing the response at the whole previous range 
of time. 

In [1], a new summation approach in technique (SAT) to find 
the transient response which does not require evaluation of all of 
these multiple reflections and leads to exact answers has been 
presented. In this paper, this technique is applied to study the 
transient response of a dipole antenna positioned as a load on 
an initially charged lossless transmission line operating as pulse 
generator. 

In Section II, a review of SAT is presented. In Section II], The 
transient input current and voltage responses for a pulse gener- 
ator loaded with a dipole antenna using SAT are investigated. 
Overall conclusions are presented in Section IV. 


II. A REVIEW OF SAT 


The summation approach in technique (SAT) to find the tran- 
sient response [1] is based on the following readily demon- 
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strated character of the Laplace transform solution for a given 
circuit configuration. Typically, this transform of v(t), V(s) is 
expanded as a sum of terms with an apparent singularity at s 
= 0. Since the limit as # — oo of v(t) = 0, the Laplace trans- 
form is, in fact, regular at s = 0. The terms of the expansion 
may contain other apparent singular points of individual terms 
$1, $2,°°-, when the transform is, in fact, an analytic function 
at these points. The net residue of the inversion integral at such 
apparent singular points of the individual terms of V(s) must 
vanish. The same may be said for z(¢) and its transform J(s). 

This character of V(s) and J(s) may be used to find expres- 
sions for the exact transient u(t) and z(t) responses consisting 
of relatively few terms (no matter what the value of time) com- 
pared to number of terms required in the usual computation by 
straightforward addition of multiple reflections. In Fig. 1, a base 
band-band pulse may be generated by discharging an initially 
charged lossless transmission line through a load. The pulse is 
initiated by closing a switch. In this system, the line is termi- 
nated in aload Z;,. Consider an initially charged section of loss- 
less transmission line, v(z,0) = Vo,0 < z < d. One way to 
get the transient response for this system is to solve the wave 
equation. The Laplace transform of the wave equation with re- 
spect to the time variable 


dV (z, 8) 


7 ea LoCo[s?V(z, 8) — su(z,0) — v’(z,0)]. 


Given the initial values v(z,0) = Vo, v’(z,0) = 0 and the 
boundary conditions I(d, s) =0, V(0, s) = —Zz,(s)/(0, s), the 
normalized transient voltage response at the load may be ex- 
pressed in the form 


1 1 
s[Zz(s) + 1] [1 — p(s) exp(-2s)] 


V(0,5)/Vo= = - 


exp(—2s) 


1 
~ s[Zz(s) +] [1 — p(s) exp(—2s)] 


(2) 


where p(s) = [Z,(s) — 1]/[Z,(s) + 1] is the reflection coef- 
ficient, Z,(s) is normalized to the characteristic impedance of 
the transmission line Zp = \/ Lo/Co, and the time variable is 
normalized so that the time for a single transit of the section of 
transmission line d LoCo = 1. 

Assume that V(0, s}/Vo is an analytic function at s = 0, and 
$1, $2,°-+, the zeros of [Z,(s) + 1]. Therefore the residues at 
each of these points in the s plane of 


Res[V (0, s)/Vo] exp(st)|s=0,s1,52,-- = 0. (3) 
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The reader is referred to [1] for the steps leading to the nor- 
malized transient voltage response for a general lumped com- 
plex load in the system shown in Fig. 1, in the form 


v(0,t)/Vo 


= S > Res 
{ [o(s) exp(—2s)]"~* exp(st) } 
s[Z1(s) + 1][1 — p(s) exp(—2s)] J 50, 51,52,-- 


+ S > Res 

{2 exp(—2s)]"~? exp(—2s) exp(st) \ 

s[Z1(s) + 1][1 — p(s)exp(—2s)] J ,20,51,82,-- 

_ pot | te(s)exp(=2s)}"™] 

L s[Zz(s)+1] | 

_ pot [exe(=2s)tels) | 

L s[Zr(s) + 1] 

_ pot | fels)exp(=28)}"™* ] 
s[Z,(s)+ 1] 


(4) 


where in the summations of residue terms n = n(t) = [(¢+1)/2+ 
1/2], n > 2, (n is the largest integer number less than n con- 
tained in the number within the brackets) and 2(n —1) < t < 
2n. In (4), L—1, means the inverse Laplace transform. 

Equation (4) is applicable for a general lumped circuit pa- 
rameter load, i.e., Zz (s) and p(s), is not function of exp(+s). 
Consequently, the transient response obtained by using (4) has 
jump points every two units of normalized values of time. If 
Z1(s) = flexp(+s)], this may lead to jump points depending 
on the nature of this function. As will be shown in the next sec- 
tion, the SAT to find the transient response for this kind of load 
can nevertheless be implemented by appropriately carrying out 
the general steps which led from (2)-(4) in the instance of a 
lumped parameter load. 


Ill. TRANSIENT RESPONSE OF A DIPOLE ANTENNA USING SAT 


In this section, SAT is applied to find the transient current and 
voltage responses of a dipole antenna modeled as a distributed 
parameter load excited by the pulse generator shown in Fig. 1, 
where the load Zy, is equivalent to the input impedance of the 
antenna. The results obtained using this technique are compared 
with those obtained by superposition of multiple reflections. 

For the system shown in Fig. 1, the normalized transient cur- 
rent response at the load may be expressed in the form 


1 


OMe =~ ele +I 


[1 — exp(—2s)] 
: (5) 
[1 — p(s) exp(—2s)] 


where Ip = Vo/Zo. 

The transmission line representation of an antenna has been 
extensively studied in [3], [13]-[15]. The authors treat the an- 
tenna as an open-circuit transmission line that is open-circuited 
at the end. Fig. 2 illustrates the representation of a center-fed 
dipole antenna as an open-circuit transmission line. The equiva- 
lent transmission line has a length equal to the half-length H of 


i(z,t) 


Fig. 1. Initially charged lossless transmission line. 


Fig. 2. The antenna as an open-circuit transmission line. 


the dipole antenna. The diameter of the antenna or transmission 
line conductor is 2a. An approximation of the input impedance 
of this antenna, [14], is 


Zin = Z.(av) coth yH = Zz (6) 


where y = a+ jf is the complex propagation constant for the 
line. Its imaginary part 2 = 27/. is the phase-shift constant, 
and its real part ~ is an attenuation constant introduced to ac- 
count for radiation and 4 is the wave length. 

Consider the particular configuration in which the average 
characteristic impedance, Z,(av), of the dipole antenna 
matched to the characteristic impedance of the pulse generator 
Zo and H = d/2. Then 


_ exp(#)exp(s) + 1 
2100) = cxp(uyexp(s) = 1 


(7) 


where 7 = ad and s is the Laplace transform variable. Conse- 
quently, the reflection coefficient will be 


p(s) = exp(—w) exp(—s). (8) 
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A. The Transient Current Response of this Dipole Antenna 
Substituting (7) and (8) into (5), one may get 
11 1 
OM = 9 5 I exp—¥)exp(38)] 
- (1 — exp(—~) exp(—s)][1 — exp(—2s)]. (9) 
Since (9) is an analytic function at s = 0, therefore the net 
residue of the inversion integral at s = 0 must vanish, 


Res[J(0, s)/Jo] exp(st)|;=0 = 0 (10) 
Expanding the term 1/[1 — exp(—w) exp(—3s)], (9) can be 
written in the form 
1(0, 8)/lo = ~5 ={1 — exp(-#)exp(-s) + exp(-¥) 
-exp(—3s) — exp(—27)) exp(—4s) + exp(—2%)) 
- exp(—6s) — exp(—3%)) exp(—7s) + exp(—37)) 
-exp(—9s) — exp(—47) exp(—10s) + +++ 
+ exp[—(n — 2)b] exp[—3(n — 2)s] 
— exp[-(n — 1) exp[-(3n — 5)s] 
+ exp[—(n — 1)] exp[—3(n — 1)s] 
— exp(-nv)exp[-(3n — 2)s}} = 5 


E 
$ 
1 


- {1 — exp(—7) exp(—s)] [1 — exp(—~)) exp(—3s)] 


- [exp(—) exp(—3s)]” + : = {exp(~2s) 

— exp(—7) exp(—3s) + exp(—) exp(—5s) 

— exp(—2z)) exp(—6s) + exp(—2z) exp(—8s) 
— exp(—3y) - exp(—9s) + exp(—3y)) 
-exp(—11s) — exp(—4w) exp(—12s) 

fee exp[—(n — 2)] exp[—(8n = 4)s] 

— exp[-(n — 1] exp[-3(n — 1)s] 

+ exp[—(n — 1)¥] exp[—@n — 1)s] 

- exp(-nii)exp(~3ns)} +5 * 

- [1 — exp(—#) exp(—s)] exp(—2s) 


: n 
Teper re) 
(11) 


where the term, 1/[1 — exp(—7) exp(—3s)][exp(—7) exp 
(—3s)]", is the remainder of the series. Let 
Fi(s) = —5 —{1- exp(—¥) exp(—s) + exp(-#) exp(-3s) 
— exp(—27)) exp(—4s) + exp(—2%)) exp(—6s) 
— exp(—37)) exp(—7s) + exp(—3%) exp(—9s) 
— exp(—4~) exp(—10s) + +--+ exp[—(n — 2)¢)] 
-exp[—3(n — 2)s] — exp[—(n — 1)y)] 
-exp[—(3n — 5)s] + exp[-(n — Dy] 
explain — 1)s] — exp(—ny) exp[—Gn — 2)s]} 
— 5 [1 — exp(—) exp(-s)| 
1 re 
j [1 = exp(—1/) exp(—3s)] [exp(-?) exp(—3s)] 
(12) 
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Fils) = 5 


{exp(—2s) — exp(—7) exp(—3s) 

+ exp(—) exp(—5s) — exp(—27)) exp(—6s) 
exp(—2%)) exp(—8s) — exp(—37)) exp(—9s) 

+ exp(—3y) exp(—11s) — exp(—4)) exp(—12s) 
+++ + exp[-(n — 2))] exp[—(n — 4)s] 

— exp[—(n — 1)#] exp[-3(n — 1)s] 

+ exp[—(n — 1)¢] exp[—(8n — 1)s] — exp(—nv) 
-exp(—3ns)} + =: 

- [1 — exp(—#) exp(—s)] exp(—2s) 


: n 
‘Toapeoenes ee eee 
(13) 


L dle 


One can find the inverse transform by evaluating the appro- 
priate inversion integral contributions due to the residues at s = 
0 only 


ResF; (s) exp(st)|,=0 
= —F{U(E) — exp(-¥ UG - 1) + exp(—-d UE 3) 
— exp(—2y))U(¢ — 4) + exp(—2y)) UC — 6) 
— exp(—3y))U(t — 7) + exp(—3y))UCE — 9) 
— exp(—4y)U(t — 10) + - + exp[—(n — 2)y] 
UE — 3(n — 2)] — exp[-(n — 1)y] 
Ue Gn—3)]+ex0l-(n— Dole 3¢0 — D] 
— exp(—n)U[E — (8n — 2)]} - gRes— 
1 
. (1 ~ exp(—v) exp(—s)] [1 _ exp(—*) exp(—3s)] 
- fexp(—7)) exp(—3s)]” exp(st)|,=0 (14) 


Resf(s) exp(st)|s=0 
= F{U(t— 2) — exp(- UE 3) t expo UE - 5) 
— exp(—2y))U(E — 6) + exp(—2y)U(E — 8) 
— exp(—37))U(t — 9) + exp(—3y)U(E — 11) 
— exp(—4y))U(t — 12) + +++ + exp[-(n — 2)y)] 
-UlE— Bn — 4)] — exp[-(n — DyJUl[t— 3@ — 1)] 
+ exp[—(n — 1)p]U[E — (8n — 1)] — exp(—ny) 


-U(E- 3n)} + SRes=[l — exp(—#) exp(—s)| 


1 
P29) TT exp) exp(—35)] 
- [exp(—7)) exp(—3s)]” exp(st)|s=o- (15) 


We must have 


Res[I(0, s)/Ip] exp(st)|s=0 
= Res Fi(s) exp(st)|,<0 + Res F2(s) exp(st)|,<0 = 0. 
(16) 
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Substituting (14) and (15) into (16) and replacing n by (n—2) 
for convenience yields the canonical result. 


; Res [l — exp(—7)) exp(—s)] i 


1 


oe exp(—3s)]"~? exp(st)|s=0 
iG Res sl — exp(—2) exp(—s)] 


-exp(—2s) 


1 
[1 — exp(—) exp(-3s)] 


fexp(-W)exp(-38)" expt) sa 
= -F{UO) — ew(-HU(E- 1) 


+ exp(—)U(E — 3) — exp(—2y))U(E — 4) 
+ exp(—27)U(t — 6) — exp(—3p)U(¢ — 7) 
+ exp(—3y)U(t — 9) — exp(—4))U(E — 10) 
+++++exp[—(n — 4))]U[t — 3(n — 4)] 


— exp(—t)) exp(—3s)] 


— exp[-(n — 3) ]U[t — (Bn — 11)] + exp[-(n — 3)¥5] 
Ult — 3(n — 3)] — exp-(n — 20) 


-U[t— Bn -8)}} + (U2) - exp(—-¥UE- 8) 


+ exp(—#)U(t — 5) — exp(—2))U(E — 6) 

+ exp(—2))U(t — 8) — exp(—3p)U(E — 9) 
+ exp(—3y)U(t — 11) — exp(—4y)U(E — 12) 
+++++exp[—(n — 4)yp]U[E — (8n — 10)] 

— exp[—(n — 3))]U[E — 3(n — 3)] 

+ exp[—(n — 3)h]UTE — (3n — 7)] — exp[-(n — 2))] 


-Ult — 3(n — 2)]} 

Excluding the remainder in (11) and taking the inverse of the 

Laplace transform, we can get the exact transient response for 
t < 3n as a superposition of multiple reflections in the form 


i(0, t)/Io = —3{U (4) — exp(—#)U(E — 1) 


+ exp(—)U(t — 3) — exp(—2H)U(t — 4) 


+++++exp[-(n — 3)P]U[E — 3(n — 3)] 
— exp[-(n — 2)y]U[E — Bn — 8) - 3 
sexpl-(n — 2WIULt— 3(n—2)]4 
sexp[-(n — I)IU[t - (Gn - 5] - 
-exp[—(n — 1)P]JU[é - 3(n - 1)]4 
-exp(—np)U[t — (8n — 2)])+ 5 

{U(E= 2) — exp(—U(E-3) 

+exp(—)U(t — 5) — exp(—2y))U(E — 6) 


NI 


NIF NIF NIE | 


+ exp(—3y)U(¢ — 11) -— exp(—4y)U(E - 
+ +++ — exp[—(n — 3)y]UTE — 3(n — 3)] 
+ exp[—(n — 3)p]U[E — (3n — 7)] 

— exp[-(n — 2)p]UTE — 3(n — 2)]} 

+ 5 exp[-(n — 2)yJUTE — B8n — 4)] 

— zexp[—(n — 1pJU[t — 3(n — 1)] 


+ 5 exp[ (n — 1)p]JU[E - Bn - 1)] 


(17) 


+ exp(—2)U(E — 6) — exp(—3p)U(E — 7) 
+ exp(—3y)U(t — 9) — exp(—4y)U(E — 10) 


+ exp(—2))U(t — 8) — exp(—3) UE — 9) 


12) 


(18) 


2.0 


1.6 


0.8 


i(0,t) /1p 


-0.8 


-2.0 


° ° ° 3 ° ° ° ° 3 
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Fig. 3. The transient current response for a pulse generator loaded with a 
dipole antenna using the infinite succession of reflections ( ) and SAT(). 


where ¢ is normalized to Ty = d\/LoCo, the travel time of the 
generator line section. 

Comparing (17) and (18), we can get the transient response 
fort < 3n in the form 


i(0,t)/Io = 5 exp[-(n — 1))]JU[E — (3n — 5)] 
— exp[—(n — 1)p]JU[E - 3(n — 1)] 
+ 5 exp(—nyp)U[E — (3n — 2)] 
+ 5exp[-(n — 1))]U[t- Bn -1)] 19) 


where n = n(t) = [(€ + 1)/3 + 1/2], (n is the largest integer 
number less than n equal to the number in the brackets),  < 
3n(t) and n(t) > 2. 

Similarly, one can get the transient current response for ¢ > 
3n in the form 


i(0,t)/Io = — 3 exp(—np)U(E — 3n) (20) 


where n = n(t) = [(£+ 1)/3 + 1/2], (m is the largest integer 
number less than n equal to the number in the brackets), ¢ > 
3n(¢) and n(t) > 0. 

We can get the transient current response using the infinite 
succession of reflections in the time domain as 
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TABLE I 
COMPARISON BETWEEN THE TRANSIENT RESPONSE OF A DIPOLE ANTENNA 
USING THE MULTIPLE REFLECTIONS CALCULATIONS AND SAT 


multiple SAT 


reflections 


: evaluation 
evaluation 


n(t) | £23n(t) | the applied the 


eq. response 


€) 
[4 [jumppoim| 6 | 2 | r<suy | a ivmppoine| « 


f2 | v<smy | ay | oo | a | 

2 20 

s | 0 -0409 | 4 
73 | -oo | wo | s | res | a» | oo | « | 
a [os | smo | am |iamp point] 1 

Ps | osu | om | oz | 1 | 


- 0.409 


19 


0 2 4 6 8 10 12 14 16 18 


t/T, 


Fig. 4. A sketch clarifies the changes in parameters necessary for use of (19) 
and (20) to find the transient response covering certain range of time. 


i(0, t)/Ip = —53{U() — exp(—p)U(t — 1) — U(E — 2) 
+ 2exp(-—p)U(E — 3) — exp(—2¢))U(¢ - 4) 
— exp(—~)U(t — 5) + 2exp(—27))U(t — 6) 
— exp(—3)U(t — 7) — exp(—2y)) 


ELS) lees (21) 

Comparison between the results obtained by using the vari- 
ation SAT represented by (19) and (20), and those obtained by 
using the infinite succession of reflections represented by (21) 
at ~ = 0.1 is shown in Fig. 3. One can see that the results are the 
same. Also, comparison between the required number of terms 
to find the transient response is shown in Table I. 

One may see the difference between the number of terms re- 
quired to evaluate the transient response at certain time using 
the infinite succession of reflections and SAT. Of course, the 
required number of terms will increase by increasing the time 
t. However, using SAT to find the transient response, the max- 
imum number of terms required as shown in Table I is 4 terms, 
whatever the increased value of time. Also in Table I we show 
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Fig. 5. The transient voltage response for a pulse generator loaded with a 


dipole antenna using the infinite succession of reflections ( ) and SAT(x). 


the alternative use of (19) and (20) according to the relation be- 
tween ¢ and n(t). A simple sketch, which clarifies the changes 
in parameters necessary for use of these equations to find the 
transient response covering certain range of time, is shown in 
Fig. 4. Of course, this sketch can be extended to cover any de- 
sired range of time. However, one can notice that the use of the 
new SAT can not cover the range of time 0 < t < 3 because 
(19) does not apply and neither does (20). 


B. The Transient Voltage Response of this Dipole Antenna 


For the same system shown in Fig. 1, starting with (2) and 
following the same steps of the variation SAT one may get an 
expression for the transient voltage response in the time domain 
in the form 


v(0,t)/Vo = 5 exp[-(n — 1)$]U[t — Bn — 5)] 
+ $exp(—ny)U[t — Bn — 2)] 


— sexp[-(n — 1LjyJU[t- Bn-1)] (22) 


where n = n(t) = [(€ + 1)/3 + 1/2], (m is the largest integer 
number less than n equal to the number in the brackets), ¢ < 
3n(£) and n(t) > 2, and 


v(0, t)/Vo = 5 exp(—ny)U(E — 3n) (23) 
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where n = n(t) = [(¢ + 1)/3 + 1/2], (n is the largest integer 
number less than n equal to the number in the brackets), ¢ > 
3n(t) and n(t) > 0. 

Also, one can get the transient voltage response using the 
infinite succession of reflections in the form 


v(0, t)/Vo = U(t) — 3{U(t) — exp(-) Ut - 1) 

+ U(t — 2) — exp(—2)U(t — 4) 

+ exp(—p)U(t — 5) — exp(—3y)U(E — 7) 
+ exp(—2y)U(t — 8) +++}. (24) 


Once more, we can see in Fig. 5 that the results at 7) = 0.1 ob- 
tained by using the variation SAT and those obtained by using 
the infinite succession of reflections are the same. However, the 
maximum number of terms, using SAT, required to find the tran- 
sient response in this case, as shown in (22), is just only three 
terms. 


IV. CONCLUSION 


The transient current and voltage responses for a pulse gener- 
ator loaded with a dipole antenna modeled as a distributed con- 
stant load using the generalized variation SAT are investigated. 
The results obtained using the SAT compared with those ob- 
tained by using the infinite succession of reflections show that 
the SAT leads to an exact answer consisting of comparatively 
few terms. The significant value of the variation SAT becomes 
evident in finding the transient response at larger values of time 
where there is no longer the need to know all multiple reflec- 
tions involved in producing the response for the whole previous 
range of time. 
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